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I. INTRODUCTION 

The study of the nuclear equation of state (EOS) and 
the behavior of nuclear matter under extreme conditions 
is crucial to our understanding of many nuclear and as- 
trophysical phenomena. With the increasing availability 
of radioactive ion-beams |1| the study of structure and re- 
actions of exotic nuclei are now possible, thus providing 
information concerning the isospin dependence of asym- 
metric nuclear matter 0, Q . Recently, the study of the 
symmetry energy for such systems has been an active 
area of interest [l, 0, H, @ . 

Most studies of the EOS involve infinite or semi-infinite 
nuclear matter and examine the dependence of the EOS 
on the parametrizations of the effective interaction as 
well as its relation to the macroscopic and macroscopic- 
microscopic models of nuclear matter Q ■ For finite nuclei 
the EOS near the equilibrium density can be investigated 
via collective observables such as the isoscalar monopole 
vibrations (breathing mode) 0]. fn addition, there are 
empirical methods relating the compressibility of finite 
nuclei to that of infinite nuclear matter Q. However, 
the behavior of the EOS for finite nuclei far from the 
equilibrium density is poorly known. 

In this manuscript we introduce a new method for cal- 
culating the zero temperature EOS and related quanti- 
ties for finite nuclei within the mean-field description of 
nuclear properties. In Section fll] we outline the general 
formalism for our calculations. Section IIIII discusses the 
application of the formalism to a sample set of nuclei 
and the results obtained. The paper is concluded with a 
summary in Section HVl 



II. FORMALISM 

In order to study finite nuclei away from their ground 
state equilibrium we take advantage of the density con- 
strained Hartree-Fock (DCHF) method 0, EJ. The 
density constraint is a novel numerical method that was 
developed in the mid 1980's and was used to provide a 
microscopic description of the formation of shape res- 
onances in light systems [ll|. Recently, we have used 
the same method to calculate heavy-ion interaction po- 
tentials from the TDHF time-evolution of nuclear colli- 
sions [12]. In the traditional constrained Hartree-Fock 
(CHF) notation, density constraint corresponds to the 



replacement 

XQ — > Xp . (1) 

The numerical procedure for implementing this con- 
straint and the method for steering the solution towards 
a specified /Oo( r ) is discussed in Refs. [13, EH- The con- 
vergence property is as good if not better than the tra- 
ditional CHF calculations with a constraint on a single 
collective degree of freedom. 

In practice, we have obtained accurate densities for 
spherical nuclei using a radial Hartree-Fock program. 
This density was then fitted to a parametric function of 
the form 

_ a (l + «ir) q 4 ,„s 

P{r > 1 + e (r-a 2 )/a 3 + 1 + e (r-a 5 )/a 6 ' W 

where ao, . . . ,a@ denote the parameters to be fitted to 
reproduce a particular density profile. In all cases the 
resulting non-linear fits were indistinguishable from the 
fitted density. Subsequently, the constraining density was 
obtained via a scale transformation of the above density 
profile 

p(r) — > p(sr) , (3) 

followed by a renormalization to produce the correct mass 
number for the nucleus under study. This scaling results 
in the compression of the bulk and stretching of the sur- 
face [8|]. Naturally, the minimum of all EOS curves occur 
at the unconstrained Hartree-Fock minimum correspond- 
ing to s = 1. The value of the scale s was generally lim- 
ited to the range (0.8, 1.2). In Fig. [T] we show the change 
in the density for various values of the scaling parameter 
s in the case of the Ca nucleus. As the primary in- 
teraction we have used the Skyrme SLy4 force ,13], with 
and without the Coulomb term, and including all of the 
spin-dependent terms. We have performed the calcula- 
tions using our new Hartree-Fock program discussed in 
Ref. [H. 

In dealing with a finite nucleus we have faced a con- 
ceptual problem of deciding what density value to use for 
plotting the density dependence of the EOS. Unless oth- 
erwise stated we have used the central density as the ref- 
erence density value for each value of the scaling param- 
eter s. Alternatively, one can choose the nuclear matter 
equilibrium density, 0.16 fm~ 3 , as the reference density 
from the density profile. The calculated values for in- 
compressibility is slightly dependent on the choice of the 



2 



reference density point due to the structure in the density 
profile. Our calculations show that this is about 10% or 
less and becomes negligible for heavy systems. One dis- 
advantage of using the central density as the reference 
value is that the equilibrium densities for different nuclei 
do not occur at the same value, thus making the calcula- 
tion of the symmetry energy, which is essentially taking 
the difference of two EOS curves along an isotope chain, 
erroneous. For this reason, in the calculation of the sym- 
metry energy we have used the density value 0.16 fm -3 
as the equilibrium density. 
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where the isospin asymmetry parameter is defined as a = 
(N — Z)/A. The a dependence of the symmetry energy 
is generally acknowledged to be 



,(/>) - S(p)a 2 + 0(a 4 ) 



(7) 



where the higher order terms in a are assumed to be 
small. Traditionally, the symmetry energy can also be 
expanded around the equilibrium density as 

L K 

4ymO) = £ S ym(po) + — (p - p ) + TTTT^ ~ + ' ' ' 

where we have defined the quantities L and K sym , which 
are related to the symmetry pressure and symmetry com- 
pressibility 



L = 3p 
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III. NUMERICAL STUDIES 

Using our approach we have first investigated the EOS 
for the 16 nucleus, with two forces, SLy4 and SkM* [H|], 
without the Coulomb interaction. We found very little 
difference between the two forces as expected since most 
modern Skyrme forces have similar nuclear matter in- 
compressibility values. For the incompressibility we find 
K A = 116 MeV. We have then repeated the same two 



FIG. 1: (Color online) Scaled density obtained for the 4S Ca 
nucleus for various values of the scaling parameter s. The 
inset shows the same quantity on a logarithmic scale. The 
force used was SLy4 without the Coulomb contribution. 



In order to extract the incompressibility coefficient, 
K\, we have expanded the EOS (binding energy per 
particle as a function of density) around the equilibrium 
density po using the expression 



E{p) _ E , K A 



= T + TM ( "-' , ° ) + 



(4) 



where Eq is the binding energy at the equilibrium density, 
and Ka is the incompressibility coefficient 



q 2 d 2 (E/A) 



dp 2 



(5) 



P=Po 



We found in practice that this expression provides an 
excellent fit except at the extreme values of the density. 
To fit the entire curve perfectly a small linear contribu- 
tion as well as a cubic term could be included. The results 
for the incompressibility is approximately 5%-10% higher 
if the full curve is used in the fit. Finally, to extract the 
symmetry energy we use the expression 
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FIG. 2: (Color online) EOS for the 40 Ca, 48 Ca, and 60 Ca 
nuclei as a function of density. Please see the comments in 
the manuscript regarding the density values. The force used 
was SLy4 without the Coulomb contribution. 



force study for the 40 Ca nucleus, yielding an incompress- 
ibility of Ka — 141 MeV. In order to investigate the 
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behavior of neutron rich systems we have performed cal- 
culations for 48 Ca and 60 Ca systems. In Fig. [5] we show 
the EOS for all of these Ca nuclei without the Coulomb 
force. Again, we stress that the density scale shown in 
Fig. [2] is determined by choosing the central density in 
the density profile of each nucleus. Consequently, the 
equilibrium value for the EOS is different for each nu- 
cleus. As we can see from Fig. [31 the 40 Ca and 48 Ca 
systems have a very similar EOS behavior, the primary 
difference being the energy shift due to the difference in 
the two binding energies. The calculated incompressibil- 
ity of 48 Ca is Ka = 155 MeV, only slightly higher than 
the one for 40 Ca. However, the situation for 60 Ca is sig- 
nificantly different, in addition to the large shift in the 
energy scale the incompressibility decreases to a value of 
Ka = 136.5 MeV, indicating a somewhat softer nucleus. 
The values for the incompressibility obtained here, per- 
haps with the exception of the 60 Ca system, are in general 
agreement with those given in Ref. [9(. 

We have repeated some of the calculations by includ- 
ing the Coulomb interaction as well. In general, EOS 
curves are shifted up due to the decrease in the bind- 
ing energy per nucleon. In addition, we see a small de- 
crease in the incompressibility modulus. For 40 Ca we 
find K A = 138.8 MeV. The incompressibility for 48 Ca 
decreases to Ka — 147.7 MeV, and this value drops down 
to K A = 123.5 MeV for 60 Ca. 
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FIG. 3: (Color online) The quantity S(p) for the 48 Ca and 
60 Ca nuclei as a function of density. 



We have also investigated the symmetry energy ob- 
tained from Eq. [51 which is essentially the difference be- 
tween the curves shown in Fig. [2] The results for the 
quantity S(p) defined in Eq. corresponding to nuclei 
48 Ca and 60 Ca having a values of 0.16667 and 0.33334, 
respectively, are shown in Fig. [3l As we observe again, 
the two systems behave very differently. Unlike infinite 
nuclear matter the curves do not cross at the equilibrium 
density since the binding energy per nucleon is different 
for different nuclei. The values of S around the equilib- 
rium density arc about 10-14 MeV, which is considerably 
lower than the estimated range for nuclear matter value 
of 32 MeV for the SLy4 force. Using the results of Fig. [3] 
we have extracted the quantities L and K syia for 48 Ca 
and 60 Ca. The results show a very large variation be- 
tween the two systems. For 48 Ca we find 66.6 MeV and 
-403 MeV for L and K sym , while for 60 Ca these numbers 
become 1.55 MeV and -193 MeV, respectively. 



IV. CONCLUSIONS 

We have introduced a new method for calculating the 
nuclear EOS for finite nuclei including all of the terms in 
the nuclear effective interaction. The calculated values 
agree well with known values obtained by other means, 
such as those deduced from experimental giant monopolc 
resonances. Much work has gone into understanding the 
effects of various components of the effective interaction 
which vanish or become very small in the infinite nuclear 
matter limit. With the increased availability of new ra- 
dioactive neutron and proton rich nuclei the study of EOS 
and symmetry energy along isotope chains of finite nuclei 
have become more urgent. We believe that the density 
constrained HF method is a step in this direction. 
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